Correlated electronic structure and chemical bonding of Ce pnictides and 7-Ce 
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We present calculated spectral properties and lattice parameters for cerium pnictides (CeN, 
CeP, CeAs, CeSb, CeBi) and 7-Ce, within the LDA/GGA+DMFT (local density approxima- 
tion/generalized gradient approximation + dynamical mean field theory) approach. The effective 
impurity model arising in the DMFT is solved by using the spin-polarized T-matrix fluctuation- 
exchange (SPTF) solver for CeN compound, and the Hubbard I (HI) solver for CeP, CeAs, CeSb, 
and CeBi. For all the addressed compounds the calculated spectral properties are in reasonable 
agreement with measured photoelectron spectra at high binding energies. At low binding energies 
the HI approximation does not manage to capture the Kondo-like peak observed for several of the 
Ce-pnictides. Nevertheless, the calculated lattice constants are in a good agreement with available 
experimental data, showing that the such a peak does not play a major role on the bonding proper- 
ties. Furthermore, the HI calculations are compared to a simpler treatment of the Ce 4/ electron as 
core-like in LDA/GGA for CeP, CeAs, CeSb, and CeBi, and the two approaches are found to give 
similar results. 

PACS numbers: 71.20.Eh, 71.15.Nc, 71.28.+d, 71.15.Mb 



h 



1 

O 
o 



> 



or 

o: 



X 



I. INTRODUCTION 

Cerium monopnictides, CeX, where X = (N, P, As, 
Sb, Bi) form in a rock-salt type crystal structure. In this 
structure the Ce atom is located on an foe Bravais lat- 
tice, with the pnictogen atoms entering octahedral voids 
of this lattice. Hence, the position of the Ce atoms in 
the Ce pnictides is the same as for the a and 7 phase of 
elemental cerium, which has attracted so much attention 
due to its isostructural volume collapse^—. The possibil- 
ity of tuning the 4/ wave function overlap between the 
cerium atoms is an opportunity in the Ce pnictides, since 
the pnictogen atoms all have different size, and hence 
they expand the foe lattice to different degrees. This of- 
fers a possibility to in a more detailed way investigate 
localization - delocalization transitions in Ce based ma- 
terials. 

The Ce pnictides have indeed been studied intensively 
in the past, due to that they display several unique and 
interesting properties, which will shortly be reviewed be- 
low. The lattice constants have been determined, and 
they suggest that the electronic structure and chemical 
bonding of CeN stand out as being different from the rest 
of the series^. This has previously been suggested to be 
caused by itinerant 4/ states for CeN&il, and the theo- 
retical interpretation- of the measured optical properties 
are in line with this reasoning^. The heavier Ce pnic- 
tides have lattice constants, and other physical proper- 
ties, which suggest a localized and non-bonding 4/ shell. 
The observed properties of the heavier Ce pnictides are 
in many ways unique, for instance, the magnetism and 
the crystal field splitting of the 4/ shell exhibit puzzling 
behaviour, as analyzed in detail in Refs. Iicl - li4l 

Another unique feature of the Ce pnictides is found 
in the excited state properties, as revealed by the x-ray 
photoelectron spectroscopy (XPS). Several investigations 
have been reported^— , and they reveal an XPS spec- 



trum which is rather complex as concerns the 4/ sig- 
nal. A feature between 2.9 - 3.2 eV binding energy is ob- 
served together with a feature closer to the Fermi level, 
at around 0.6 eV binding energy. The exact positions 
of these two distinct features vary from pnictide to pnic- 
tide, see e.g. Ref. [IB and [T3- The presence of these two 
peaks is at first glance surprising, having in mind that a 
localized f 1 configuration of the heavier Ce pnictides, is 
expected to result in one spectral feature corresponding 
to a f 1 — > f° transition. However, it has been pointed 
ou t 17 ' 22 that two screening channels of the final state core 
hole are possible for this excitation, a well screened and 
a poor screened channel, which results in two distinct 
spectral features. The presence of two peaks in the 4/ 
spectrum is a property not only of Ce pnictides, but also 
other Ce compounds (see, for example, Refs. [23Tj27l ). 

On the theoretical side, electronic structure calcula- 
tions have been published for CeN, analyzing the energy 
band structure and the optical properties^. Results of 
this work are based on the LDA/GGA and, essentially, 
correlation effects are taken into account only weakly. 
A recent calculation 2 ^ for CeN focused on the electronic 
structure and the lattice dynamics. In this work the 4/ 
states were considered as delocalized, and the calcula- 
tions used both the LDA/GGA level of approximation 
as well as the so called GW— approximation. For the 
heavier Ce pnictides a theory of the crystal field split- 
ting was proposed^, which successfully rep roduced the 
observations of Ref. [jjj- The work of Ref. [30 was a combi- 
nation of model considerations and first principles theory, 
in which the localization of the 4/ wave function was en- 
forced in the calculation, and then the hybridization to 
other valence states was estimated from this electronic 
configuration. 

Previous DMFT investigations of Ce pnictides are pre- 
sented in Refs. [3lTj33l . wherein the Anderson impurity 
problem was solved in a finite-U extension of the non- 
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TABLE I. Muffin- Tin radii (in atomic units) of the cerium 
and pnictogen atoms used in the present calculations. 





CeN 


CeP 


CeAs 


CeSb 


CeBi 


R M T(Ce) 


2.45 


2.75 


2.83 


2.90 


2.89 


-Rait(X) 


2.10 


2.40 


2.50 


2.90 


2.91 



TABLE II. Experimental binding energies (AE) and DC 
terms for the 4/ states of cerium in the Hubbard I approxi- 
mation. 





CeP 


CeAs 


CeSb 


CeBi 


AE (eV) 


2.95 


3.05 


3.0 


3.2 


PW91, DC (Ry) 


0.335 


0.308 


0.262 


0.247 


PBE96, DC (Ry) 


0.344 


0.316 


0.270 


0.255 


AM05, DC (Ry) 


0.356 


0.328 


0.283 


0.269 



crossing approximation (NCA), including some lowest- 
order crossing diagrams. 

In the present paper we revisit the problem of the elec- 
tronic structure of the Ce pnictides, and the localization- 
delocalization phenomena of the 4/ shell, in light of re- 
cent developments of treating electron correlations di- 
rectly in electronic structure calculations. In particular, 
we investigate how well dynamical mean field theory re- 
produces the lattice parameters and spectral properties 
of these materials, and we have used the implementation 
reported in Refs. I34T437I for these purposes, with the HI 
and SPTF solvers. 



II. COMPUTATIONAL DETAILS 

The calculations presented in this work are based on 
the density functional theory (DFT ) 38 i 39 combined with 
the dynamical mean- field theory 4 ^ method. The DFT 
calculations were done in the local-density approximation 



TABLE III. Calculated and experimental lattice parameters 
of cerium pnictides in atomic units. 





CeN CeP CeAs CeSb CeBi 


PW91, 4/-itinerant 
PBE96, 4/-itinerant 
AM05, 4/-itinerant 


9.38 10.73 11.01 11.72 11.96 
9.55 10.94 11.26 12.00 12.27 
9.42 10.78 11.09 11.81 12.07 


PW91, 4/-core 
PBE96, 4/-core 
AM05, 4/-core 


9.79 11.08 11.31 11.92 12.10 
9.98 11.32 11.59 12.24 12.57 
9.86 11.17 11.41 12.03 12.20 


PW91+DMFT 

PBE96+DMFT 

AM05+DMFT 


9.38 11.09 11.32 11.95 12.13 
9.56 11.30 11.55 12.22 12.44 
9.43 11.19 11.43 12.06 12.25 


Experiment 


9.49 11.20 11.49 12.13 12.26 



with exchange-correlation functional PW91— , and in the 
generalized gradient approximation using the PBE96 46 
parametrization, as well as with the AM05 functional 4 ^, 
within the full-potential linearised muffin-tin orbitals 
(FP-LMTO) approacbJMs^. The LDA/GGA+DMFT 
calculations were done by means of the implementation 
presented in Refs. 13414371 . with the total energies evalu- 
ated following Ref. |35|. To solve the Anderson impurity 
model arising in the DMFT part, the SPTF solver^&^i 
was used for CeN, whereas the Hubbard I solve r 36 ! 40 ' 52 
was used for CeP, CeAs, CeSb, and CeBi. 

The LDA/GGA+DMFT scheme was applied to the 
4/ state of cerium using two types of the correlated 
orbitals 3 ^: orthogonalized LMTO subspace (ORT) and 
muffin-tin only correlated subspace (MT). All these ORT 
and MT orbitals are localized, and form orthonormal 
sets. The ORT set consists of functions obtained from 
the FP-LMTO method with the LDA/GGA Hamiltonian 
for the correlated states. The MT correlated subspace is 
defined only into the muffin-tin region, and outside of 
the muffin-tin spheres the MT orbitals are taken equal 
to zero by definition. Properties of these subspaces are 
analyzed in the works published in Refs. HH and l53l . In 
the present calculations the MT set was used for CeN, 
and the ORT set for all other Ce pnictides. 

The LDA/GGA calculations included spin-orbit cou- 
pling in the scalar-relativistic approximation, and only 
within the muffin-tin spheres. The radii of the muffin-tin 
spheres for the different compounds were kept constant 
while changing the lattice parameters, and their values 
are reported in Table HI For the Brillouin zone integra- 
tion the tetrahedron method was used with a 16 x 16 x 16 
k-mesh. It corresponds to 85 k-points in the irreducible 
wedge of the Brillouin zone with the total number of k- 
points in the Brillouin zone equal to 2048. For the charge 
density and potential angular decomposition inside the 
muffin-tin spheres the value of maximum angular mo- 
mentum was taken equal to Z max = 8. Three kinetic en- 
ergy tails were used, corresponding to {—0.3, —2.3, —1.6} 
Ry. Finally all calculations were non-spin polarized, ex- 
cept CeN which was converged to a ferromagnetic solu- 
tion. Magnetic properties were in fact considered at the 
DMFT level, and the SPTF solver works better when 
starting from a magnetic solution^ 4 -. 

Calculations without DMFT in the pure 
PW91/PBE96/AM05 approximations were done in 
two ways: the cerium 4/ electrons were treated either as 
core or as itinerant states. In the latter case two kinetic 
energy tails were considered for the basis. 

For the LDA/GGA+DMFT calculations it is neces- 
sary to choose an adequate value for the Hubbard and 
(Hund) exchange parameters U and J. For all Ce pnic- 
tides except CeN we have used U=7 eV and J=0.95 eV, 
which is compatible with previous calculations 31 consid- 
ering that our method is also based on a fully rotationally 
invariant U-rnatrbs 3 ^. Conversely the 4/ states in CeN 
have basically an itinerant character, as discussed more 
in detail in the results section. This implies that the 
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screening of the Coulomb interaction is much more ef- 
fective and the value of U reduced in comparison to the 
other Ce pnictides. Following these guidelines and hav- 
ing in mind that the SPTF solver tends to overestimate 
correlation effect s 35 ' 55 we have used U=2 eV, which is 
consistent with its applicability range. Higher values of 
U have also been studied but these results have not been 
reported here as they gave the same physical picture as 
\J—2 eV. The value of the Hund's exchange is less depen- 
dent on the nature of the screening in a material, and 
therefore we have used J=0.95 eV for CeN as well. 

To calculate the energy versus volume curves in the 
Hubbard I approximation, it is important to correctly 
treat the double counting term Hue included in the lo- 
cal problem^ 5 .. The Hoc term is used to remove the 
wrong contributions of local Coulomb interaction already 
contained in the DFT Hamiltonian projected on the cor- 
related orbitals (4/ states) and used in the DMFT part. 
The double counting term was first calculated from the 
HI solver at experimental lattice constants, and then 
kept fixed for all the other values. In order to fix the 
main occupied peak at the correct position the photoe- 
mission data available from experiment s 16 ' 18 ! 19 ' 56 were 
used. The experimental binding energies and double 
counting term values are presented in Table |TT] for the 
PW91/PBE96/AM05 + HI approximations. 

CeN did not present these problems since in the SPTF 
solver the double counting correction can be directly re- 
lated to a fully hybridizing self-energy by simply remov- 
ing the orbitally averaged S(0), separately for each spin 
channel. When evaluating the local interacting Green's 
function with the SPTF solver was important to consider 
a fully renormalized propagator to obtain a conserving 
approximation in Baym-Kadanoff sensed— . 
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FIG. 1. Experimental and calculated lattice constants of Ce 
pnictides. The theoretical values for CeN are obtained with 
the SPTF solver and the AM05 DFT functional. The Hub- 
bard I solver together with the AM05 DFT functional was 
used for all other pnictides. 



with observed lattice parameters. As Fig. [5] shows, the 
difference between the two functionals, when combined 
with DMFT, is very small and both parametrizations re- 
produce observations with good accuracy. In Fig. [2] we 
also show results from more time efficient calculations, 
treating the 4/ states as localized core states and using 
the AM05 DFT functional. Moreover for CeN a calcu- 
lation with the same functional and 4/ states treated as 
itinerant is showed. 

The DFT calculations with 4/ electrons treated as core 
states for CeP, CeAs, CeSb and CeBi, give a good agree- 
ment with the experimental data. On this level of ap- 
proximation the 4/ states do not contribute to the chem- 
ical bonding. The Hubbard I calculations also treat the 
4/ states as essentially non-bonding and atomic like, al- 



III. RESULTS 

The calculated lattice parameters are presented in 
Table IIII1 where we also show available experimental 
dat a 5 ' 61 ~ ^£, A selection of these results are plotted in 
Fig. [1] In this figure we show experimental values of 
lattice constants and values obtained with AM05 DFT 
functional combined with the SPTF solver for CeN and 
with the Hubbard I solver for all other pnictides. An 
overall good agreement between experimental data and 
calculations is observed, both concerning the trend as 
well as the absolute values. 

In Fig. [5] we present the calculated lattice constants 
from different approximations. Results are given for CeN 
using the SPTF solver with AM05 and PW91 DFT func- 
tionals as well as for pure AM05 DFT calculations with 
the 4/ electron of Ce atom considered as itinerant. For all 
other pnictides, results are given for the Hubbard I solver 
with AM05 and PW91 DFT functionals. The choice of 
showing in Fig. [5] parametrizations from the AM05 and 
PW91 functional was made since these two approxima- 
tions give the best and, respectively, the worst agreement 
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FIG. 2. Calculated lattice constants of Ce pnictides. The 
theoretical values for CeN are obtained with the SPTF solver 
and two DFT parametrizations, AM05 and PW91. The Hub- 
bard I solver with both the AM05 and PW91 DFT functionals 
was used for all other pnictides. These data are compared to 
calculations using localized 4/ (core) states with the AM05 
and PW91 functional. For CeN the plain AM05 DFT results 
with 4/ valence electrons are also reported, emphasizing how 
such a description is much better than the assumption of lo- 
calized core states for this material. 



FIG. 3. (Color online) Experimental and calculated spectrum 
of CeN. The blue points label the XPS data and the red points 
mark the BIS spectrum from Ref. [68l . where the relative in- 
tensities of the XPS and BIS spectra have been tentatively 
adjusted. The factor ( x 10) comes directly from Ref. Ir38l . and 
we have kept it to make a better comparison between the ex- 
perimental XPS data and the theoretical curve. The solid line 
denotes the total spectral function using the SPTF solver for 
the Ce 4/ states. 



beit a small degree of hybridization is present in such 
kind of calculations. All data in Fig. [T]and Fig. [5] point 
to that the 4/ states do not contribute to the chemical 
bonding in the heavier Ce pnictides. 

In contrast, the CeN compound is a system which 
seems best described having delocalized 4/ states as it 
follows from the comparison of our calculated values and 
the experimental data. The DMFT calculation with the 
SPTF solver and the pure AM05 DFT calculation with 
itinerant 4/ states reproduce the observed lattice con- 
stant with good accuracy. Conversely if one tries to as- 
sume a localized 4/ electron in the form of a core state the 
results become more in line with the other Ce pnictides, 
but the agreement with experiments worsen. These facts 
point to that the 4/ states have an important role in the 
chemical bonding, which results in a significantly smaller 
lattice constant of CeN in comparison to the heavier Ce 
pnictides (FigfT|). The delocalized nature of the 4/ elec- 
tron makes it not suitable for the Hubbard I approxima- 
tion or any other approach in the strong coupling limit, 
unless one wants to apply semi-empirical workarounds. 
For example, in a previous LDA+DMFT study^L based 
on the NCA solver, the 4/ hybridization function in CeN 
had to be rescaled of a factor 2 to avoid an unphysical 
occupation number of the 4/ shell. 

To give further evidence to this fact we observe that 
the calculated spectral properties of CeN are in rather 
good agreement with experimental data, when consider- 
ing the 4/ states as essentially itinerant with some degree 
of correlation resulting from the SPTF solver. This com- 
parison is presented in Fig. [31 where both occupied and 
unoccupied states from the calculations are compared to 
experimental XPS and BIS data. Note that the calcula- 
tion of the occupied states reproduce the spectral feature 
between 3 and 5 eV binding energy, stemming mostly 
from N p-states, whereas the feature at the Fermi level is 



FIG. 4. (Color online) Experimental and calculated spec- 
tral properties of the 4/ states for CeP. The XPS data from 
Ref.El are marked as blue circles. The red circles label the 
resonant spectrum of CeP, while the red dashed line labels the 
extracted photoemission spectrum of 4/ states from Ref. 
The black solid line is the calculated spectral function of the 
Ce 4/ states with the Hubbard I solver. 



dominated by the Ce 4/ states, which form a band-like 
feature pinned at the Fermi level. This results in that the 
unoccupied states have a 4/ feature at the Fermi level, 
which is different from the situation for the heavier Ce 
pnictides, as discussed below. The unoccupied states of 
CeN also have a broader feature at 5 - 6 eV energy, which 
is dominated by Ce bd states which hybridize with ligand 
p states. 

In Fig. [4] we present the CeP 4/ spectral function, 
which sometimes is loosely referred to as the correlated 
density of states (DOS) to relate it to the band-structure 
theory. In this figure we make a c omp arison to the ex- 
perimental data from Refs. [HI and [Tj- As one can see, 
the Hubbard I approximation describes the deeply local- 
ized 4/ feature in a proper way, although due to life time 
broadening and finite spectrometer resolution the exper- 
imental feature is much wider compared to theory. Note 
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FIG. 5. (Color online) Total spectral function calculated with 
the Hubbard I solver (black solid line), compared with exper- 
imental data from Ref. [l7l . We also present the projections 
for the Ce 5d and 4/ states, and the P 3p states (dashed 
lines). The vertical arrow labels the pd peak which is close to 
the experimental non resonant photoemission data taken at 
hu = 112 eV (red solid line). 
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FIG. 6. (Color online) Experimental and calculated spectral 
properties using Hubbard I approximation of CeAs. The total 
spectral function is labelled by the black solid line, the Ce 
4/ projection is marked by the black dashed line, while the 
red dashed line marks the experimental resonant 4/ emission 
from Ref. IvH . The red solid line represents the experimental 
data from nonresonant emission (formed by As 4p and Ce 5d 
states). 



that the Hubbard I approximation is not able to capture 
the well screened feature at the Fermi level, and indeed 
the experimental feature can not be reproduced with the 
present calculations. 

The experimental results for the non resonant fre- 
quency of 4/ Ce states are compared with the total and 
partial theoretical spectral functions in Fig. [SJ For this 
photon energy the main contribution to the observed 
peak is from cerium 5d states and phosphorus 3p states^. 
The experimental peak-width and position are repro- 
duced by the Ce 5d and P 3p partial spectral functions, 
as Fig. [5] shows. 

A similar picture occurs for the other heavy pnictides, 
for which we in Figs. [51 El [3] compare calculated spectral 
features with experimental data taken on and off reso- 
nance. As discussed, the 4/ feature close to the Fermi 
energy (the well screened peak) is not reproduced in our 
Hubbard I calculation. However, going from CeP to CeBi 
the experimental 4/ peak at the Fermi level decreases 
in intensity and becomes almost negligible in the CeBi 
case. For this pnictidc the Hubbard I approximation 
hence gives the best description of the spectral features. 

In Figs. El HI [71 and [5] we also show the unoccupied 
states. It is clear that the 4/ dominated spectral fea- 
tures are due to an / — > f 2 transition, in which multi- 
plet effects are clearly visible. Since we have not found 
experimental BIS data for the heavy Ce pnictides our 
theory must be viewed as a prediction. We note that a 
calculation based on the LDA or LDA+U approximation 
would not result in spectral features displaying multiplet 
effects, and an experimental verification or rebuttal of 
the BIS data shown in Figs. [5j [BJ El [S] would give valu- 
able insight into which approximation is best for these 
systems. On general grounds however, one would expect 
for these materials that the Hubbard I approximation is 
better than LDA or LDA+U. 

A similarity can be observed between the heavier Ce 




E - E F ,(eV) 

FIG. 7. (Color online) Total spectral function of CeSb calcu- 
lated with the Hubbard I solver (black solid line), compared 
with the experimental data from Ref. [l6l . We also present 
projections for the Ce 4/ and 5d states, and the Sb 5p states 
(dashed lines). The pd-peak is close to the experimental data 
from nonresonant photoemission (red circles). 



pnictides and 7-Ce. In the photoelectron spectrum 
of Ce both the 7-phase and the a-phase display the 
well screened 4/ feature and the poor screened 4/ fea- 
ture, which for 7-Ce is located at about 2.0 eV binding 
energ y 18 ' 69 . 7-Ce is a localized 4/ system, and there- 
fore a Hubbard I calculation for 7-Ce should reproduce 
the lattice constant with equal accuracy as it does for 
the Ce pnictides. Our calculations for 7-Ce actually con- 
firm this expectation and we obtain an equilibrium lat- 
tice constant of 9.56 a.u., 9.95 a.u, 9.78 a.u., respectively 
for PW91+HI, PBE96+HI, and AM05+HI. These values 
should be compared to the experimental lattice parame- 
ter which is equal to 9.75 a.ui^il, showing again that the 
AM05+HI approach results in the best agreement with 
experiment. For 7-Ce we have also calculated the bulk 
modulus to investigate if the lack of the well-screened 
4/ feature in the spectrum leads to a deficiency in more 
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FIG. 8. (Color online) Experimental and calculated spectral 
properties of CeBi using the Hubbard I approximation. The 
total spectral function is labelled by the black solid line, while 
the Ce 4/ contribution is marked by the black dashed line. 
The red dashed line and the red solid line respectively mark 
the experimental resonant 4/ emission from Ref. Q3 and the 
nonresonant emission, formed by the Bi 6p and Ce 5d states. 
The blue dashed line labels the resonant 4/ emission from 
Ref. GJ. 
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sensitive ground state properties than the lattice con- 
stant. We obtain a bulk modulus of 37 GPa, 31 GPa, 
and 27 GPa, respectively for PW91+HI, PBE96+HI, and 
AM05+HI. These values are in line with the GGA+U 
result of 34 GPa£^, and overestimate the experimental 
bulk modulus, which is equal to 19 GPaZi. A better 
theory based on a more sophisticated approach would 
probably be able to correct this disagreement with the 
experiments. 



IV. CONCLUSIONS 

We conclude that the LDA/GGA+DMFT approach 
used in the present work with the SPTF and Hubbard I 
solvers gives very good values of lattice constants and 
acceptable spectral properties of Ce pnictides. The Hub- 
bard I solver is not able to capture the experimental 
Kondo-like feature at the Fermi level observed for CeP, 
CeAs and CeSb, while for CeBi this peak is almost ab- 
sent, and only a small shoulder is visible. This defi- 
ciency is inner to the method used, and may potentially 
be remedied with more sophisticated techniques such as 
quantum Monte- Carlo. However, in this case, technical 
problems of different nature may be introduced and re- 
lated to the statistical noise of the Green's function. For 
reasonable computational resources the occupation num- 
bers and the total energies would be less stable, and the 
spectral function would suffer of the uncertainty related 
to the maximum entropy method for analytical continua- 
tion. Finally it would also be problematic to include the 
spin-orbit coupling without arbitrarily discarding parts 
of the hybridization function. 

It is interesting, nevertheless, that despite the inability 
of the Hubbard I approximation (or the SPTF solver) to 



capture the well-screened Kondo-like peak observed for 
some Ce pnictides the lattice constants are reproduced 
with a very good accuracy. This implies that the mecha- 
nism behind this feature docs not contribute much to the 
chemical bonding. More sensitive properties as, e.g., the 
bulk modulus would probably show a signature of this 
discrepancy between theory and experiment, as verified 
in this paper for 7-Ce. It would be of a great interest to 
approach this problem with more refined methods, and 
this will be the aim of our future research. 

Finally we also note that the lattice constants for the 
heavier Ce pnictides, in which the 4/ states are essen- 
tially localized and non-bonding, can be well reproduced 
from a theory where these states are treated as core- 
states where they are not allowed to hop from lattice 
site to lattice site, and where hybridization to all other 
states is absent. Hence, a good alternative to treat the 4/ 
states of localized 4/ electron systems is to simply con- 
sider them as core like, which is an approximation that 
offers great computational speed with decent accuracy. 
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